Exchange instabilities in electron systems: Bloch versus Stoner Ferromagnetism 
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We show that 2D and 3D electron systems with the long-range Coulomb electron-electron inter- 
action could develop ferromagnetic instabilities due to strong exchange effects at low densities. The 
critical densities in both 2D and 3D systems at which the magnetic instability, which could either be 
of Stoner type (second-order) or of Bloch type (first-order), are higher than the dispersion instability 
critical density where effective mass at the Fermi surface diverges. We discuss the theoretical as well 
as experimental implications of the ferromagnetic instability at low electron densities, particularly 
in low-disorder semiconductor-based two-dimensional systems. 
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I. INTRODUCTION 

With the development of high mobility semiconduc- 
tor based 2D electron systems such as Si MOSFETs 
and GaAs HIGFETs, where extremely low carrier den- 
sities and very high quality can be achieved, a number 
of recent experiments has been carried out to measure 
fundamental physical quantities such as magnetic sus- 
ceptibility 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 : 10 ! 11 ! 12 ! 13 and effective massii in 
low density 2D electron systems at low temperatures 
(< lOOmK). These experiments have triggered many 
theoretical as well as experimental studies on 2D and 3D 
electron systems. One important purpose of these stud- 
ies is to obtain the ground state phase diagram for such 
systems. It is well known that the zero temperature 2D 
and 3D electron systems can be characterized by a sin- 
gle dimensionless interaction parameter r s = E e - e / Ek 
(where E e ^ e denotes the interaction energy and Ek the 
kinetic energy) , which depends only on the density (n) of 
the system with r s cx n -1 / 2 (n -1 / 3 ) in 2D (3D) systems. 
For small r s , electron systems can be well described by 
Landau's Fermi liquid theory. (Actually, there is an ex- 
ception even to this as at exceptionally low temperatures 
an interacting Fermi liquid undergoes a Khon-Luttinger 
superconducting transition, which we ignore for our pur- 
pose since it is of no physical relevance.) In the limit of 
large r s , the system tends to reduce its interaction en- 
ergy E e _ e at the cost of higher kinetic energy by forming 
into an electron crystal (the so called Wigner crystal^), 
which has been established by Monte Carlo studies in 
both 2Di& (where the crystallization transition happens 
at r s ~ 32 — 42) and 3D systemsii (where the transi- 
tion happens at r s ~ 55 — 75). How the system behaves 
in between the above mentioned two limits of r s is not 
yet clear. It has been widely accepted, and also sug- 
gested by Monte Carlo studies, that there may exist a 
ferromagnetic phase (of either fully or partially polarized 
spins) in the intermediate r s region. Many other theoret- 
ical studies for such high r s value region typically start 
with a more or less arbitrary assumption of a particular 
ground state symmetry of the system. Among the vari- 
ous model of exotic interaction-driven electronic ground 



states, charge or spin density wave states, various super- 
conducting states, glassy or clustered ground states have 
been much discussed in the literature. 

From a theoretical perspective, it is important to study 
the evolution of the Fermi liquid state starting from the 
weakly interacting (small r s ) regime. As r s increases, 
all the single particle properties of the system are in- 
creasingly renormalized by interaction effects. The ques- 
tion is whether the Coulomb interaction renormalization 
brings about certain instabilities eventually at some large 
r s , and hence changes the ground state symmetry. The 
existence of degrees of freedom related to spin and mo- 
mentum makes it natural to consider the possibility of 
a magnetic instability or an instability that is related 
to the dispersion of a quasiparticle, which we call "dis- 
persion instability" for short. In spite of a great deal of 
past theoretical work investigating the magnetic instabil- 
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sion instability possibilities 3 ^ 3 ^ in both 2D and 3D elec- 
tron systems using a large variety of different theoreti- 
cal techniques, a unified theoretical treatment of the fer- 
romagnetic instability, starting from the weak-coupling 
Fermi liquid ground state with interaction effects intro- 
duced systematically, is still lacking. This absence of a 
unified picture has left many unanswered important ques- 
tions. First, what is the order of any possible ferromag- 
netic instability? It could be a second-order transition 
(the so-called Stoner instability) caused by the continu- 
ous divergence of the susceptibility as the critical density 
or r s -value is approached from the weak-coupling side. 
Or it could be a first-order transition (the so-called Bloch 
ferromagnetism) which could happen abruptly at a spe- 
cific r s -value. Both have been predicted and studied in 
the literature, but their inter-relationship has not been 
clarified. Second, what is the relationship between mag- 
netic instability and dispersion instability? Whether the 
possible divergence of the interacting spin susceptibility 
X* , which can be written as \* — 9* m *i where m*(g*) 
are the effective mass (^-factor) of the system, is caused 
by the divergence of m* (dispersion instability) or g* , or 
both? In order to answer these questions, we provide 
in this paper a comprehensive picture of the Bloch ferro- 
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magnetic instability, Stoner ferromagnetic instability and 
the dispersion instability of the 2D and 3D electron sys- 
tem using the method of Random Phase Approximation 
(RPAp.. 

The structure of the paper is the following: in Sec- 
tion [HI we discuss briefly the background for Bloch and 
Stoner ferromagnetism; in lllll we describe our theoretical 
formalism; in IIVI we present and discuss our theoretical 
results for ferromagnetic instability in interacting 2D and 
3D quantum Coulomb systems; we conclude in [V] with a 
discussion of experimental implications and related open 
questions. 



II. BACKGROUND 

The possibility of a density driven ferromagnetic tran- 
sition in an interacting electron system was first sug- 
gested by Bloch- 5 ' more than 75 years ago. Bloch's ba- 
sic idea, essentially a Hartree-Fock mean field theory, 
remains fundamentally valid even today. The idea is 
that at high density the electron system would be para- 
magnetic in order to optimize the kinetic energy cost 
(which is high in a high-density quantum fermionic sys- 
tem) whereas at low density the system should sponta- 
neously spin-polarize itself into a ferromagnetic ground 
state in order to optimize the exchange energy arising 
from the Pauli principle and Coulomb interaction. For 
an electron gas (in a positive jellium background) it is a 
straightforward exercise to write down the total Hartree- 
Fock energy per-particle as a sum of the non-interacting 
kinetic energy and the (Fock) exchange energy due to 
unscreened Coulomb interaction at T = as 
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for 2D systems, where n = rij + ni is the total num- 
ber density of electrons; £ = (n^ — n^/n is their spin- 
polarization (or magnetism) density; a g n = ( ^ 7rr^ ) — 1 
and a 2 B n = (irr 2 ) -1 define the dimensionless interac- 
tion parameters r s in 2D and 3D respectively (with 
(ib = ft 2 /(me 2 ), the Bohr radius); and the energy is 
measured in Rydberg units (i.e. e 2 /as)- It is easy to 
see that the above Hartree-Fock energy expressions lead 
to a first-order ferromagnetic transition (the "Block fer- 
romagnetism" ) at r s = tb where tb — 5.45 (3D) and 
tb — 2 (2D), i.e. E(( = 1) ferromagnetic state is lower 
(higher) in energy than E(( = 0) paramagnetic state for 
r s > (<)tb- We refer to such energy-difference-driven 
abrupt (first order) transition as Bloch ferromagnetism 
in the rest of this paper. 



The Stoner ferromagnetic instability-- refers to the di- 
vergence of the spin susceptibility and hence a second- 
order continuous magnetic phase transition from a para- 
magnetic (r s < rst) weak-coupling side to a ferromag- 
netic side (r s > rst)- The simplest model to con- 
sider is, following Stoner's original work, a zero-range 
delta-function like interaction of strength l F (a constant 
in momentum space) between the electrons, leading to 
an interacting static long-wavelength spin susceptibility 
(in the dynamical Hartree-Fock approximation) given by 
X* lx = (1 — DqI)^ 1 , where D = D(Ep) is the electronic 
density of states at the Fermi energy. This immediately 
leads to the Stoner criterion for ferromagnetic instability 
defined by a divergent x* lx when 1 — DqI = 0. Since 
the 2D density of states is a density-independent con- 
stant, this instability criterion does not lead to a mean- 
ingful condition in 2D unless we arbitrarily define / to 
be the Coulomb interaction strength at Fermi wavevec- 
tor, whence the Stoner instability criterion leads to un- 
physically low r s values for the ferromagnetic instability 
given by rst — 0.7 (2D) and rst — 1-5 (3D), which are 
absurdly small revalues and are unrealistic. Of course, 
in real electron systems the electron-electron interaction 
is the long-range Coulomb interaction, and therefore the 
simple Stoner instability criterion, defined by DqI = 1 
where / is an effective short-range interaction strength, 
is inapplicable. But the basic idea of the ferromagnetic 
Stoner instability, defined by a continuous divergence of 
X*( r s) as r s is increased, still applied. We refer to the 
ferromagnetic transition defined or characterized by a di- 
vergence of the interacting susceptibility as the Stoner 
instability. 

In real electron liquids, the exchange-only Hartree- 
Fock approximation considered above for the Bloch in- 
stability is inadequate because correlation effects (i.e. en- 
ergy contributions beyond Hartree-Fock) are known to be 
extremely important, and must be included in the ener- 
getic considerations. Similarly, the interacting suscepti- 
bility must be calculated for the real Coulomb interaction 
in the system, not for a hypothetical zero-range interac- 
tion, in order to obtain a better estimate of the Stoner 
instability criterion. In the rest of this paper we consider 
both Bloch and Stoner ferromagnetic instabilities using 
better many-body approximations for Coulomb electron 
liquids, namely RPA. 



III. THEORY 



In this work we follow the notation of Refs. 36 3 7ll4ll 
Within RPA 36 i 37 i 38 i 41 i 42 i 43 i 44 , the Coulomb contribution 
to the ground state energy of a jellium electron system 
with long-range Coulomb interaction can be denoted by 
the Feynman diagrams shown in Fig. ^V. The quasi- 
particle energy is then obtained by E^ — SEa/Sn^, 
where nk is the distribution function at momentum k. 
The second order derivative of the total ground state 
energy is referred to as Landau's interaction function: 
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FIG. 1: The RPA Feynman diagram for: (A) the Coulomb 
interaction contribution to the ground state energy; (B) Lan- 
dau's interaction function; (C) Dynamically screened inter- 
action. The circles are polarization bubbles, the thin wiggly 
lines are the bare Coulomb interaction, and the solid lines the 
noninteracting electron Green's function. 



/(k, k') = 5 2 Eq I Sn^Sn^' , represented by the Feynman 
diagram shown in Fig. ^3. Graphically, taking the nk 
variational derivative of a quantity simply means that 
one cuts one solid line of the Feynman diagram and takes 
the external momentum and frequency to be on-shell (i.e. 
lu = k 2 /2m — Ep with m the band electron mass and 
Ep the Fermi energy). We emphasize that the RPA as 
shown in Fig. ^ necessarily implies that this on-shell self- 
energy approximation is used for calculating the quasi- 
particle energy dispersion Ek and Landau's interaction 
function /(k, k') since all energy and momenta in Fig.Q] 
correspond to the noninteracting system. Thus RPA self- 
energy approximation necessarily implies an on-shell ap- 
proximation (Fig. ^V) as emphasized by Rice— a long 
time ago. Solving the full Dyson equation using the off- 
shell renormalized energy, as is sometimes done in the 
literature, would be completely inconsistent within RPA 
as it would incorrectly mix various perturbative orders. 

Following Hubbard's notation^, the ground state en- 
ergy per particle Eq/N with N the particle number 
can be written as E G /N = E K /N - w(0)/2 + E c /N, 
where Ek is the kinetic energy, v(0) — J v q d d q/ {2ir) d 
is the interaction energy at zero separation with the 
bare Coulomb interaction v q — 2ire 2 / q for 2D and v q — 
4ire 2 /q 2 for 3D, and Ec is the Coulomb contribution to 
the ground state energy (both exchange and correlation) 
which can be denoted as in Fig.QjV. Note that the singu- 
larities in v(0)/2 and Eq/N cancel out with each other. 
In the 2D system, it is easy to show^ that RPA leads to 
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xdx I du[\ne(x, u) — e(x, u)], (3) 



where e(q,u>) is the dynamical dielectric function, g s is 
the spin degeneracy (g s — 2 for paramagnetic states and 
g s = 1 for fully spin-polarized state), E ex is the exchange 
part of the ground state energy. Note that in Eq. © 
we subtracted a term e(q,tj) from ln[e(<7, a>)] in order to 



handle the ultraviolet divergence in the integration. Sim- 
ilarly for 3D we have 
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x 2 dx I du[lne(x, u) — e(x, u)\. (4) 



It is convenient to convert all the expressions in terms 
of the dimensionless units r s . The definition of r s is as 
r s = 1/ '(akpdB) where Uf is the Fermi momentum and 
as = l/(me 2 ) is the Bohr radius. In 2D a = v / gJ/2 and 
in 3D a = (2g B /9ir) 1 ^ 3 . Also we choose h = 1 through- 
out, which makes wavevector and momentum (as well as 
energy and frequency) equivalent. In these notations, it 
is easy to show that E ex /N = —8ar s / (2>n)Ep for 2D, and 
E ex /N = —3ar s /(2ir)Ep for 3D. In the actual calcula- 
tion, the integration in Eqs. © and (0J can be performed 
on either real or imaginary axis. By examining the r s and 
g s dependence of the ground state energy, we study the 
Bloch magnetic instability of the electron system. The 
integrals in Eqs. © and (0} are the correlation contri- 
butions not included in our Hartree-Fock considerations 
of Bloch ferromagnetic in Sec. [n] Note that g s = 1(2) 
corresponds to ( = 1(0) in Sec.lTTl 
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FIG. 2: (a) The ladder-bubble series for the interacting sus- 
ceptibility with the bold straight line the interacting Green's 
function and the bold wavy line the dynamically screened 
interaction; (b) the noninteracting susceptibility; (c) the 
Dyson's equation for the interacting Green's function in terms 
of the noninteracting Green's function and the self-energy; (d) 
the self-energy in the leading-order expansion in the dynam- 
ical screening; (e) the Dyson's equation for the dynamically 
screened interaction in terms of the bare Coulomb interaction 
(thin wavy lines) and the polarization bubble; (f) a charge 
fluctuation diagram which does not contribute to spin sus- 
ceptibility. 

We investigate the Stoner instability by calculating the 
magnetic susceptibility \* within RPA which is repre- 
sented by the Feynman diagram showed in Fig. [21 Direct 
calculation of these diagrams turns out to be difficult for 
the long ranged Coulomb interaction. However, at T = 0, 
Landau showed that x* can be equivalently expressed as 
the following equation^: 



m 
m* 



fe(0)do, 



(5) 
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where x is the Pauli spin susceptibility, f e {9) = /e(k, k') 
with k and k' on-shell; k 2 /2m = k' /2m(= Ep) is the 
exchange Landau's interaction function, 6 is the angle 
between k and k', do the element of solid angle along k' 
times cos in 3D and dO in 2D, and vj = l/(27r) 2 in 2D 
and w = kp/(2ir) 3 in 3D. Similarly, the Landau theory 
expression for the effective mass m* is^ 



to 

= 1 — w 

to* 



f(9)do. 



(6) 



In Eq. JSJ, f e (9) is the exchange part of the Landau's 
interaction function, which is represented by Fig.^3(a). 
This spin independent term is responsible for the differ- 
ence between the ratio x/x* an d to/to*—. An equivalent, 
and easier way to derive the effective mass is through cal- 
culating quasiparticle self-energy and obtaining its mo- 
mentum derivative. The self-energy within RPA can be 
written as2& 



£(k,w) 



d d q f dv v q 



-G (q+k,i/+w), (7) 



(2ir) d J 2m e(q, v) 
where d = 2 or 3 is the dimension of the system, and 



G (k,w) = 



1 - nA£k) 
i ~ £k + iO^ 



(8) 



is the bare Green's function, where np is the Fermi dis- 
tribution function, £t = k 2 /(2m) — Ep. It is shown^i 
that the integration along real axis in the expression of 
self-energy (Eq. J7J)) can be deformed onto imaginary 
axis, which avoids the singularities along the real axis 
and makes the integration easier. The contour deforma- 
tion also breaks the expression of self-energy into sepa- 
rate terms that correspond respectively to contributions 
from the spin-dependent and spin-independent part of 
the Landau's interaction function shown in Fig. ^ and 
is very useful for us to derive the expression for suscep- 
tibility as shown below. The expression of the real part 
of the self-energy can then be written as 

Re £(k, w) = - J ^v q Q(2mw + k 2 F - |q - k| 2 ) 



d d q 



w Re 



(2ny Uq - e (q,£ q _ k - W ) 
0(2muj + k 2 F -\q- k| 2 ) - Q(k 2 F - |q - k| 2 ) 



d d q 
(2^p 



dv 
— i 
2tt 



e(q,iu) 



1 



W - Cq+k 



(9) 



The effective mass is derived from the expression of the 
real part of the quasiparticle self-energy by m/m* = 
1 + {m/k F )^-Re E(k, £k)|ft=fci^- Combining this with 
Eq. JSJ, we have 



/ f e (0)do 



m d 
kp dk 



ReE(k,a)U= 



(10) 



It is not difficult to show that the second term of 
Eq. 10 accounts for the contribution from the spin- 
independent exchange Landau's interaction function 
/ e (k, k') (Fig.np(a)), and therefore the term w J f e (9)do 
in Eq. I|10|l exactly cancels the momentum derivative of 
the second term in the self-energy Eq. ©. Hence the 
expression of x/x* onr y contains contributions from the 
k derivatives of the first and third term in Eq. JjJJ • Af- 
ter converting all the expressions in terms of the dimen- 
sionless parameter r s , and using 2kp, AEp, 2m as the 
momentum, energy, and mass units, we obtain the ex- 
pression for the 2D magnetic susceptibility as 
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(11) 



where A — x — x 2 — 
larly for 3D we have 



i,B = 2xu, R = VA 2 + B 2 . Si 
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dx 



du 



e(x, iu) 



x [ln(F/G) - 2C/F + 2D/G] 



(12) 



where C = 1 - q, D = 1 + q, F = C 2 + u 2 , G = D 2 + u 2 . 
Note that in the expressions of the dielectric function 
e(x, u), x — qj (2kp) and u~ujj (AEp). 



IV. RESULTS 

In Fig. we present the calculated magnetic suscep- 
tibility as a function of r s for both 2D and 3D systems, 
together with the calculated g-factor g* — x* l m * an d 
the effective mass to* . It is clear from Fig. [21 that both 
2D and 3D systems experience Stoner ferromagnetic in- 
stabilities, characterized by the divergence of magnetic 
susceptibility as the density decreases. It is important to 
note that this Stoner instability (i.e. divergence of %*) 
does not arise from an effective mass divergence since the 
to* divergence happens at much lower densities. In other 
words, g* and m* both diverge, with the divergence of g* 
occurring at a lower r s value. For the 2D system, x* an d 
g* diverges at r s ~ 7.3 while to* diverges at r s ~ 18.1. 
For the 3D system, x* an d 9* diverge at r s ~ 18.7 while 
to* diverges at r s ~ 49.9. 

We emphasize that both to* and g* divergences actu- 
ally happen independently and are completely unrelated 
phenomena. On the other hand, it is also worth mention- 
ing that the g* divergence does have some quantitative 
effect on the to* divergence. After g* diverges, the sys- 
tem becomes a ferromagnetic liquid, and the momentum 
distribution is different and the Fermi energy increases. 
This change results in a small increase in the critical r s 
value where to* diverges. In fact for g s = 2 paramagnetic 
systems, to* diverges at 16.1 for 2D and 47.8 for 3D2&, 
in contrast to g s = 1 case where to* diverges at 18.1 for 
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FIG. 3: (A) and (C): Calculated renormalized spin suscepti- 
bility x* IXi effective mass m*/m and g-factor g*/g. For 2D 
system, \* an d 9* diverge at r 3 ~ 7.3 while m* diverges at 
r s ~ 18.1. For 3D system, \* an d 9* diverges at r s ~ 18.7 
while m* diverges at r s ~ 49.9. (B) and (D): inverse suscep- 
tibility shows x* diverges at r s ~ 7.3 for 2D and 18.7 for 3D 
systems. Note that the x* an d g* are calculated for paramag- 
netic systems (g a — 2) while m* are for ferromagnetic systems 

(5a = 1). 



valley degeneracy = 1 



(A) 
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FIG. 4: Quasi-2D effects on magnetic susceptibility diver- 
gence. (A): GaAs quantum well system. (B): Si-inversion 
layer. 



2D and 49.9 for 3D. But this effect is rather small and is 
of no particular significance. 

In real experimental systems, the value of r s at which 
X* diverges should be influenced by many factors (even 
within our RPA many-body approximation scheme). 
Here we consider the valley degeneracy g v and the finite 
width effect on the magnetic susceptibility divergence in 
semiconductor-based realistic 2D electron systems. The 
effect of g v is exactly the same as the effect of g s on the 
system, and therefore can be easily incorporated. For 
the finite width effect, we introduce a form factor to the 
Coulomb interaction, followi ng the standard procedure 
described in detail in Refs. I45l46l Using appropriate 
semiconductor parameters, we obtain the susceptibility 
in GaAs quantum wells and Si-inversion layers, plotted 
in Fig. It is clear from Fig. 0] that multi- valley degen- 



eracy and finite width both suppress the divergence of 
the susceptibility renormalization, and make the critical 
r s value for x* divergence considerably larger than the 
strict 2D results. 
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FIG. 5: (A) and (C): the RPA ground state energy per particle 
for polarized (g s = 1) unpolarized (g s = 2) cases as a function 
of r s for 2D and 3D systems, calculated using Eq. @ and 
Eq. (B) and (D): details around the magnetic instability 
at r 3 ~ 5.5 for 2D and 17.8 for 3D systems. 

The magnetic susceptibility is a thermodynamic Fermi- 
surface property. As mentioned before, another way 
of studying the magnetic instability is to compare the 
ground state energy of the system for polarized and un- 
polarized states at different electron densities. This is the 
Bloch ferromagnetism discussed in Sec. [HI Our results 
(Fig. |3J) of RPA ground state energy for fully polarized 
and non-polarized electron states (using Eqs. © and Q) 
in both 2D and 3D electron systems show similar char- 
acteristic. When r s is very small (or electron densities 
high), both systems prefer non-polarized paramagnetic 
states. As r s increases to a certain critical value (r s ~ 5.5 
for 2D and r s ~ 17.8 for 3D) the ground state energy for 
the fully polarized electron state actually becomes lower 
than the non-polarized states. This clearly indicates that 
the system undergoes a Bloch type ferromagnetic insta- 
bility due to the Coulomb electron-electron interaction in 
a low density 2D electron systems. Note that the critical 
r s for the Bloch instability is substantially higher in the 
RPA theory (increasing from 2 to 5.5 in 2D and 5.45 to 
17.8 in 3D) than in the Hartree Fock theory due to the 
inclusion of correlation energy. 

Comparing Figs. and [S] we conclude that, at least 
within our well-defined RPA ring-diagram many-body 
approximation scheme, the sequence of instabilities (as 
density decreases) the theory predicts for both 2D and 
3D electron liquids is the following: Bloch ferromag- 
netism (i.e. an abrupt first order magnetic transition) 
at r s = 5.5 (2D); 17.8 (3D); Stoner ferromagnetism char- 
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acterized by a continuous divergence (i.e. a second order 
magnetic transition) of the interacting g-factor and of the 
susceptibility at r s — 7.3 (2D), 18.7 (3D); the dispersion 
instability associated with the continuous divergence of 
the quasiparticle effective mass at r s — 18.1 (2D), 49.9 
(3D). Of course, comparing first and second order tran- 
sitions is not particularly meaningful since their origins 
are fundamentally different, and a first order transition 
may always preempt a second order transition as seems 
to happen in Coulomb electron liquids with Bloch ferro- 
magnetism always happening (both in 2D and 3D) at a 
slightly lower r s values (5.5 versus 7.3 in 2D, and 17.8 
versus 18.7 in 3D) although the difference in the critical 
r s values for the two transitions (less than 10% in 3D 
and about 25% in 2D) is sufficiently small so that both 
Bloch and Stoner ferromagnetism remain of experimental 
interest. 

We also note that, by definition, x* = g*m* (and 
X = gm) , and therefore the divergence of the interacting 
susceptibility could be caused either by a divergence g- 
factor or a diverging effective mass. This issue has been 
much discussed and debatediSiiiSiiiiiSii^ in the recent 
experimental literature on 2D semiconductor-based elec- 
tron systems, where low-density divergence of both 2D 
susceptibility and effective mass has been reported. All 
we can say is that our theoretical results are only consis- 
tent with the susceptibility divergence as arising from the 
divergence of the interacting g- factor, not the effective 
mass, since the g-factor divergence occurs at much lower 
r s -values, r s ~ 7.3 (18.7) in 2D (3D) for g* j g divergence 
versus r s ~ 18.1 (49.9) in 2D (3D) for m* /m divergence. 
We add that in realistic quasi-2D semiconductor system 
(our Fig. 0J) the susceptibility (as well as effective mass) 
divergence occurs at substantially higher r s values due 
to the considerable softening of the Coulomb interaction 
from its strict 2D form due to the finite- width effect. 



V. DISCUSSION 

Both Bloch and Stoner instabilities imply that 2D 
and 3D electron systems interacting via the long range 
Coulomb interaction undergo a T = ferromagnetic 
quantum phase transition from a high-density paramag- 
netic state to a low-density ferromagnetic state either 
through a first order (Bloch) transition or a second-order 
(Stoner) transition (with a continuous divergent suscep- 
tibility) at a critical r s -value. Given that the critical r s 
value(s) for the transition(s) we obtain within our RPA 
many-body theory is rather large (i.e. r s 3> 1), we do 
not expect our predicted r s parameter for ferromagnetic 
transitions in 2D and 3D electron systems to be reliable. 
But the basic trends, such as the sequence of instabilities 
(i.e. Bloch followed by Stoner followed by the dispersion 
instability with diverging mass as the density is lowered) 
or the suppression of the transition to much lower densi- 
ties in quasi-2D systems, should be valid, in general. In- 
deed quantum phase transitions predicted by RPA (or for 



that matter even by the simpler Hartree-Fock approxima- 
tion) are always found to occur in the numerical quantum 
Monte Carlo (QMC) simulation albeit at higher r s values. 
This is certainly true for electron liquid ferromagnetic 
instabilities. For 3D system, QMC simulation predicts 
that the system undergoes a possibly second-order phase 
transition and electron spins become partially polarized 
at r s ~ 60 by Ref. or r s ~ 15 - 25 by Ref. [H and 
as r . in crease to r s ~ 100 by Ref. \^ or r s ~ 35 — 45 by 
Ref. [13, the system becomes a fully polarized ferromag- 
netic state. For 2D case, QMC simulation predicts^ that 
as density approaches r s ~ 25, the system undergoes a 
first-order transition into a fully spin-polarized ferromag- 
netic state. Theoretically it is possible to obtain the spin 
susceptibility through calculating the second derivative 
of the ground state energy with respect to £ at C = 0. 
However, in reality too much error is introduced when ob- 
taining the susceptibility this way in QMC simulations, 
and therefore it is difficult to ascertain the order of the 
ferromagnetic transition in QMC numerical simulations 
- our RPA theory predicts the first-order Bloch instabil- 
ity to occur first as the density is being lowered. Another 
thing to be noticed about QMC simulation results is that 
the ferromagnetic transition r s value predicted by differ- 
ent groups differs by a large factor from each other (see, 
for example, Ref. and Ref. l47f) . which shows the large 
amount of error introduced by such simulations due to 
the small energy difference between spin polarized states 
and spin unpolarized states (the two density dependent 
energy curves are almost parallel when they cross each 
other) and different choices of trial wave functions. 

Much has been written about the validity of the RPA 
many-body approximation at low carrier densities (r s > 
1). We have little to add to this issue beyond the de- 
tailed discussion we already provided in our recent pub- 
lications^Siai. We want to emphasize that, although RPA 
is exact in the r s — > limit, it is by no means a theory 
based on an r s expansion - it is a self-consistent field the- 
ory based on an expansion in the dynamically screened 
interaction which should be qualitatively valid for all r s 
below the Wigner crystallization of the electron liquid. 
In fact, RPA is found to be quantitatively valid in 3D 
systems^ at metallic densities (r s ~ 3 — 6) and in 2D 
systems for r s up to 10 — 15 where comparison with ex- 
periment has been carried ouli^. Often the error in the 
calculation arising from other effects (e.g. finite temper- 
atureSi, finite quasi-2D width, band structure, etc.) turn 
out to be larger than that included in the RPA approx- 
imation, and therefore improvement beyond RPA (short 
of a full-fledged QMC calculation) becomes meaningless. 
One can try to "improve" upon RPA by including local 
field corrections to the dynamical electron polarizability 
(i.e. bare bubble of RPA) which, in some crude manner, 
simulates the incorporation of higher-order vertex cor- 
rections in the theory. But such local field corrections 
are uncontrolled, and probably inconsistent, since many 
diagrams in the same order are typically left out. We are 
therefore unconvinced that the inclusion of local field cor- 
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rections in the theory is necessarily an improvement on 
RPA. The great conceptual advantage of RPA is that it 
is a well-defined approximation that is both highly phys- 
ically motivated (i.e. dynamical screening) and theoreti- 
cally exact in the high-density (r s — + 0) limit. Attempted 
improvement upon RPA through the arbitrary inclusion 
of local field correction may neither be theoretically jus- 
tifiable nor more reliable. Keeping these caveats in mind 
we have carried out our exchange instability calculations 
using the Hubbard local field corrections^ 3 -, and we find 
no qualitative changes from the RPA results presented 
in this paper. The critical r s values for the occurrence 
of the ferromagnetic instability change somewhat in the 
presence of the local field correction, but this is a re- 
sult without any significance since the precise values of 
critical r s are expected to be not particularly accurate 
in any of these theories. The fact that the basic quali- 
tative conclusions about the various instabilities do not 
change in the presence of local field corrections demon- 
strates the qualitative robustness of our RPA-based re- 
sults. Another point to note is that the fact that RPA 
predicts the existence of divergence of certain physical 
quantities at certain critical densities does not necessar- 
ily imply that RPA or Fermi liquid theory fails at that 
density. For example, RPA predicts the divergence and 
then negative values for the compressibility 38 at densi- 
ties higher than ferromagnetic transition densities, and 
negative compressibility has indeed been observed in ex- 
periments^ on 2D electron systems as predicted by RPA 
calculations. Of course the total compressibility of a sys- 
tem cannot be negative, but just the electronic part of 
the compressibility can be negative as predicted by RPA 
and as is routinely observed in 2D electron systems for 
r s > 3. There can be no doubt that if 3D electron systems 
with large enough r s values (r s > 6 according to RPA 38 ) 
are found they would routinely have negative electronic 
compressibility as well! It is certainly true that RPA be- 
comes a progressively poorer approximation as density 
decreases and perhaps detailed QMC calculations should 
be carried out to test the validity of our RPA-based pre- 
dictions presented in this paper. 

In discussing possible experimental implications of our 
results, we note the great recent experimental inter- 
est in the literature on the possibility of a density- 
driven ferromagnetic transition in semiconductor-based 
2D carrier systems where very low carrier densities 
(r s < 10) can be achieved in rather high-quality sam- 
pleo 1 ! 2 . 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 . There are recent experi- 
mental claimsiSi&i^ of the observation of a low density 
susceptibility divergence in Si MOSFETs for r s ~ 7 — 10. 
Although the experimental susceptibility behavior as a 
function of density (or, r s ) looks similar to our theoret- 
ical results in Figs. and 0] we are skeptical about the 
significance of this agreement. There are several reasons 
for our skepticism. First, the experimental claimed sus- 
ceptibility divergence occurs at far too high a density 
(r s ~ 7) compared with the theory where we find the 
RPA susceptibility divergence in realistic Si MOSFETs 



(Fig. to be occurring at r s w 23. This RPA predic- 
tions for critical r s (~ 23) is most likely the lower bound 
- any real susceptibility divergence is expected to occur 
at higher r s values (r s > 23). Second, the experimen- 
tal divergence of x* Ix = g*m* /gm has been claimed 1 * 2 * 3 - 
to be arising from an effective mass divergence, not the 
g-factor divergence as we find in our theory. Our the- 
oretical effective mass divergence, in fact, occurs at a 
critical r s more than twice as large as the corresponding 
X* critical r s . In fact, our quasi-2D effective mass diver- 
gence 3 ^! occurs for r s > 40! Third, there has been no 
experimental evidence for the existence of a low-density 
ferromagnetic phase such as hysteresis, remanence, etc. 
If there is indeed a ferromagnetic transition, one should 
be able to observe ferromagnetic behavior at densities be- 
low the ferromagnetic critical density (i.e. for r s values 
larger than the point of x* divergence). No such direct 
ferromagnetic behavior has ever been observed in a low- 
density 2D system casting serious doubts on the claims 
of the observation of a 2D ferromagnetic transition. A 
very recent extremely careful and detailed measurement 
of 2D susceptibility in a high-mobility n-GaAs system^ 
finds no divergence in x* Ix up to r s sa 12, calling into 
question the earlier claims of susceptibility divergence in 
Si MOSFETs at lower r s values. In addition, a direct 
thermodynamic measurement of 2D susceptibility in Si 
MOSFETs also does not find a ferromagnetic instabil- 
ity. What is clear is the observed strong enhancement of 
X* Ix (and m*/m) as a function of increasing r s which is 
consistent with our theoretical findings. But the actual 
existence of a low-density electron liquid ferromagnetic 
transition has not be established experimentally in our 
opinion. 

Finally we discuss some of the earlier literature that is 
of relevance to our work. Our calculation of the ground 
state energy for polarized and unpolarized states par- 
tially confirms the numerical results of Rajagopal et. 
alm^, who also considered the possibility of partially 
polarized states, and found that for certain range of 
electron densities, the system prefers a partially polar- 
ized ferromagnetic state in the 3D system. Similar re- 
sults were also derived using the Quantum Monte Carlo 
metho d 17 ! 47 . In this paper we did not present our calcu- 
lation results (which confirms the results of Rcf. 52 ) for 
the Bloch instability associated with the partial spin- 
polarization because this will not help our understand- 
ings of the relation between Stoner, Bloch, and disper- 
sion instabilities, which is the main purpose of this pa- 
per. In 3D, partially spin polarized states is preferred 
energetically for the density region of r s between ~ 14 
and ~ 18, which is right before the fully polarized fer- 
romagnetic region as the density is decreased^. This 
suggests that our understanding of the relation between 
the three kinds of instabilities will not be affected by the 
consideration of partially spin-polarized states. In 2D 
systems, partial spin-polarization does not occur^ 7 ***** 2 -, 
and thus this issue does not arise at all for our 2D cal- 
culations, which is the main focus of our work. For the 
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magnetic susceptibility, there have been earlier RPA cal- 
culations in 2D 29 i 30 i 31 i 53 i 54 and 3D2£, and QMC calcula- 
tions in 2D^£. Only ShastrySS predicted a susceptibility 
divergence in 3D systems, and our results confirm his con- 
clusion. No previous work considered the inter relations, 
among Bloch instability, Stoner instability and dispersion 



instability. Our work is the only work in the literature 
connecting all these density-driven electron liquid insta- 
bilities within one coherent theoretical framework. 

This work is supported by NSF, DARPA, ONR, LPS, 
and ARO. 
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